o 

(N 

C 
3 



- 1—1 

x 



An eccentrically perturbed Tonks- Girardeau gas 

J. Goold 

Centre for Quantum Technologies, National University of Singapore, Singapore, 
Singapore 



■ M. Krych, Z. Idziaszek 

Faculty of Physics, University of Warsaw, 00-681 Warsaw, Poland 

OJO 1 T. Fogarty, Th. Busch 

Department of Physics, University College Cork, Cork, Ireland 

a 

1 Abstract. We investigate the static and dynamic properties of a Tonks- 

^j-i Girardeau gas in a harmonic trap with an eccentric 5-perturbation of variable 

strength. For this we first find the analytic eigensolution of the single particle 
problem and use this solution to calculate the spatial density and energy profiles 
of the many particle g ELS 3iS cL function of the strength and position of the 
perturbation. We find that the crystal nature of the Tonks state is reflected 
in both the lowest occupation number and momentum distribution of the gas. 
' O ' As a novel application of our model, we study the time evolution of the the 

Ch ' spatial density after a sudden removal of the perturbation. The dynamics exhibits 

O , collapses and revivals of the original density distribution which occur in units of 

O i the trap frequency. This is reminiscent of the Talbot effect from classical optics. 
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1. Introduction 

One of the standout achievements of ultracold atomic physics in the past decade 
has been the realisation of controllable and clean environments for the simulation of 
quantum manybody systems [JJ. Cold atoms are ideal systems for such simulations 
and offer many desirable features, two of which are particularly unique. The first of 
these is the ability to tune the scattering length and hence the interaction in quantum 
gases. This is achieved by means of Fcshbach resonances [3J. The second is the ability 
to change the dimensionality of the system and in particular to generate periodic 
geometries for atoms by means of optical lattice potentials [3] . Both of these separate 
developments, or their combination, allows experimentalists to enter regimes of strong 
correlation and explore complex condensed matter models. For instance, using optical 
lattice potentials the superfluid-Mott insulator phase transition has been observed by 
changing the lattice depth [3j[4] . In addition, by loading a BEC into a two-dimensional 
lattice and utilising strong confinement in the transverse direction, the achievement 
of lower dimensional systems was made possible and new quantum phases can now 
be explored with an incredible degree of precision. In the case of two dimensions, 
the famous Bcrczinskii-Kosterlitz-Thouless phase transition was observed recently [5] 
and for ultracold gases in one dimension one of the crowning achievements was the 
realisation of a unique model from many body theory [B] - the strongly correlated 
Tonks-Girardeau gas of hard core bosons [5]- 

Since the first creation of a Tonks-Girardcau gas there have been further 
experimental studies focusing on the absences of thcrmalisation due to the intcgrability 
of the underlying Hamiltonian [9] and more recently on the density dynamics in the 
presence of an outcoupled impurity [10]. This second work has opened the door 
to studying the transport properties of impurities in ultracold quantum gases and 
therefore highlights the importance to understand the fundamental properties of a 
Tonks-Girardeau gas in the presence of an impurity. A number of studies in this 
direction already exist in the literature, in which the impurity is modelled using 
pseudo-potential approximation [Til El [JJ2 [HI [HI HE] . These studies have primarily 
focused on fixed position perturbations and in this work we introduce a versatile 
analytical model which can be used to describe the Tonks-Girardcau gas in the 
presence of a perturbation of arbitrary strength at any position in a harmonic trap. In 
addition to describing a static impurity, this model can be interpreted as the limiting 
case of an split, asymmetric double well trap which may be realized using a sharply 
focused laser beam which is detuned from the atomic transition. 

The paper is organised as follows: in Sec. [2Jwe briefly remind the reader of the 
Tonks-Girardeau gas and its solution by the Fermi-Bose mapping theorem. We then 
outline the eigensolution of a harmonically trapped single particle in the presence 
of an off-centre (5-perturbation in Sec. [3] and show in Sec. [4] how this can be used in 
combination with known techniques to characterise the groundstate of a many-particle 
gas. We then study the dynamics of the single particle density after a sudden removal 
of the perturbation as well as for a timc-of-flight measurement. Finally, we conclude. 

2. The Tonks-Girardeau gas 

Optical lattices and atom chips allow to create trapping potentials that are tight 
enough in the transversal direction to freeze out all dynamics in these degrees of 
freedom [17] . A gas of N bosons in such a potential can then be approximated by the 
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onc-dimcnsional Hamiltonian 

+ 9id^2S{\x 1 - Xj \) , (1) 

i<j 

where m is the mass of the particles, V ext is the trapping potential and gn> is a ID 
coupling constant which is derived from the rcnormalisation of the three-dimensional 
scattering process gm = 4 f n °^ p (a± — Cg^d) -1 HE]- Here a± is the trap width and 
C is a constant of value 1.4603.... This Hamiltonian describes an inhomogeneous 
Lieb-Liniger gas, which in the strongly repulsive limit, gm — > oo, can be solved 
by using a mapping to an ideal and spinless fermionic system [6]. This procedure 
is known as the Fermi-Bose mapping theorem and it can be used to show that 
the local density and correlation functions of this strongly correlated system are 
equivalent to the corresponding quantities of a non-interacting spin polarized Fermi 
gas. This strange equivalence is due solely to the dimensionality of the system - as the 
repulsive interactions become stronger, the particles are no longer free to overlap, thus 
mimicking the Pauli-exclusion principle in configuration space. The essential idea is 
that one can then treat the interaction term in eq. ((TJ) by replacing it with a boundary 
condition on the allowed bosonic wave-function 

ty B (xi,X2, . . . ,x n ) = if \xi-Xj\=0, (2) 

for i ^ j and 1 < i < j < N. This is simply the hard core constraint which says that 
no probability exists for two particles ever to be at the same point in space. Such a 
constraint is automatically fullfillcd by calculating the wave-function using a Slater 
determinant 

I (JV-l.JV) 

ty F (xi,X2,...,x N ) = -t=- det ip n (xj) , (3) 
V N\ (« j')=(o,i) 

where the ip n are the single particle eigenstates of the ideal system. This, however, 
leads to a fermionic rather than bosonic symmetry, which can be corrected by a 
multiplication with the appropriate unit antisymmetric function [S] 

A = sgn(x l -xj) , (4) 

l<i<j<N 

to give — A$f. Knowing the single particle eigenstates is therefore a pre-requisite 
to being able to apply the mapping mechanism. 



K 2 d z 

2m dx"i 



+ V ex t(x n ) 



3. Single particle problem 

The number of models in many particle quantum physics for which analytical solutions 
exist is extremely limited and often restricted to one spatial dimension |19| . Here we 
will describe one such model by studying the Hamiltonian ([1]) with 

Vext (x) = \^mu?x 2 + j5(x - d) , (5) 

and giD — > oo, which describes a Tonks- Girardeau gas in a harmonic trap of frequency 
uj which is split asymmetrically by a tunable (5-perturbation of strength 7 a distance d 
from the centre of the trap. The limit d = for this model is well known and has been 
studied extensively in recent years |TTJ [13l HH [20] . In order to use the prescription 
of the mapping theorem one needs to first solve for the single particle eigenfunctions 
to build the Slater determinant © . In this section we will study the eigcnsolution of 



An eccentrically perturbed Tonks- Girardeau gas 4 

the single particle problem in detail. In fact, this problem is equivalent to the relative 
problem of two interacting atoms in seperate harmonic traps which has recently been 
numerically solved in the three dimensional case . We show below that this problem 
is analytical in one dimension. 

The single particle Schrodingcr equation for the potential (0 is given by 

'-ft 2 d 2 1 \ 

- — 3-5- + -mu 2 x 2 + 7<5(x - d) ip„(x) = e v ^ v {x) , (6) 
2m dx 1 2 J 

where the energies are given by e u = (v + ^)haj. Let us rescale all of the units in terms 

of the undisturbed (7 = 0) trap length l^ = J and energy tuj. In this way the 
equation (|6]) can be rewritten as 



dx 2 +V+ 2 ~T ~ k6 ( x ~ d ) ) M x ) = 0, ( 7 ) 

where k = jIq/(Hlj) is the re-normalised strength of the <5-barrier. On either side 
of the 5-function the solutions of the differential equation (J7]) are parabolic cylinder 
functions D l/ (x), which vanish for x — > 00, but diverge for x — > —00. We can therefore 
write the solution piecewise as 

ip u (x) = ij}i jV (x)6(d -x) + tp r>u {x)e{x - d) , (8) 

with 

iprA x ) = N+D v {x) and ^i, v {x) = N_D v {-x) . (9) 

and 9(x) being the Heaviside function. The condition of continuity of these solutions 
at the position of the <5-function 

N + D u {-d) = N„D v {d) , (10) 

together with the solution of the Schrodingcr equation 

ijj v {x)dx + / V(x)ip u (x)dx = e u / ip„{x)dx , (11) 

J d—e J d—e 

with 

V(x) = -x 2 + kS(x - d) , (12) 

leads to a transcendental equation which determines the energy eigenvalues as a 
function of both, k and d 

N+D'J-d) + N_D' u {d) - 2nN + D u {-d) = . (13) 

The derivatives of the parabolic cylinder functions can be calculated using the 
recurrence relation [52] 

Dl(z)=uD v _ 1 (z)-^zD I/ (z), (14) 

and for D v (d) ^ we find from eq. (HUJ) the relation iV_ = N + D v {-d)/D v (d). 
Substituting these into eq. (fLlj) we therefore find 

D v _ x {-d)D v {d) + D v _ 1 (d)D u (-d)^ = 2KD v {-d)D u {d) . (15) 

This equation determines the eigenenergies for the solutions that are nonzero at the 
position of the delta potential. To find the ones for which D v (d) = we can see from 
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Figure 1. The energy spectrum of harmonically trapped particle in the presence 
of 5-like perturbation at position d = 1 and d = 4, for different strengths, re, of 
the perturbation. 




Figure 2. The energy spectrum of harmonically trapped particle in the presence 
of (5-like perturbation with strength re = 1 and re = 10, for various positions of the 
perturbation. 



cq. |T0]) that D v (-d) = 0, provided that N + ± and iV_ ^ 0. Then eq. ((T3j) reduces 
to 

N + vD v _x{-d) + N-vD v ^(d) = 0. (16) 

where the derivatives have been determined using eq. (|14[) . The above equation 
together with the normalization of the total wave function determines N + and 7V_ 
for the solutions that vanish at the position of the delta potential. These solutions 
are independent of k, and they correspond to the harmonic oscillator wave functions 
that vanish at x = d. 
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Figure 3. The lowest four cigenfunctions of a particle in a harmonic potential 
with a repulsive 5-function located at d = 0.5 of strength re = 10 



In Figs.Q]and[2]we show the energy spectrum as a function of d and k, respectively. 
The presence of the perturbation introduces a non trivial structure in the spectrum and 
in general leads to an increase of the state's energy. The four lowest lying eigenfunction 
for d = 0.5 and k = 10 are shown in Fig. [3] It can be seen that the presence of delta 
potential results in an abrupt change in the slope of the wavefunction at the position 
of the delta function and that the states are already highly localised on one side of 
the impurity for k = 10. 

4. Static and Dynamic Properties 

4--1- Single Particle Density Profile 

The single particle density is one of the most important observables for ultracold 
quantum gases. In the Tonks regime one can obtain it, even time-dependently, from 
the spectrum of underlying single particle Hamiltonian as |23j 



where we have adopted the convention of labeling the first N eigenfunctions as 
n = 0,1,2,... N. 

In Fig. [4] we show this single particle density for a gas of 20 particles in a trap with 
an impurity of strength k = 10 at three different positions in the trap. As expected 
the ^-impurity creates a significant local dip in the density and has only minimal effect 
at larger distances. The enhanced oscillations which are present around the position 
of the impurity are analogous to the famous Friedel oscillations which occur around 
an impurity in the surface charge density of a homogeneous electron gas [24] ■ 




N-l 




(17) 
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Figure 4. The single particle density of a harmonically trapped Tonks-Girardeau 
gas with 20 particles in equilibrium with a barrier of strength k = 10, located at 
positions d = 0, 3, —6. 

Let us in the following investigate the dynamical properties of such a gas by 
examining a non-equilibrium situation created by sudden removal of the impurity. In 
order to compute the time-dependent density of cq. (|17|) one needs to employ a time 
dependent basis. We obtain this basis numerically using the split operator method in 
the unperturbed harmonic trap. Alternatively one may simply employ the well known 
propagator for the harmonic oscillator to get the time dependent basis |25| . This can 
then be used to calculate the single particle density of eq. (fTT)) . 

In Fig. [5] we show p(x, t) following the sudden removal of an impurity of strength 
K = 10 located at d = —6. One can see that the density dip formed by the impurity 
vanishes almost instantaneously, however a mirror image of the dip appears after half 
a trap period, t = tt/w, and then again dissapears followed by a complete revival after 
one trap period, t = 2tt/uj. This effect is analogous to the Talbot effect from classical 
optics where periodic rcfocusing of a diffraction grating is expected to occur in the 
near field of a transmitted wave. In this situation the (5-function represents the most 
trivial form of a diffraction grating. A natural question to ask is why is such an effect 
occurring in a strongly correlated many body system as the Talbot effect is a coherent 
single particle effect? The simple answer can be found by recalling that the system 
can be mapped onto free fermions, for which the single particle density is simply the 
sum of the squares of the single particle eigenfunctions, with each one undergoing its 
own coherent unitary evolution. In this picture, the occurrences can be explained by 
noting that all N eigenfunctions superimpose in phase again after every trap period 
and one finds that the density profile at odd multiples of tt/lj is a mirror image of the 
initial density profile at t = 0. It is also worth noting that in between this revivals 
the density shows an interesting fine structure. This is shown in a close up of the 
density for the time period < t < tt in Fig. [6] When the impurity is removed, the 
mattcrwave readjusts to the profile of a harmonically trapped gas by a relaxation of 
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Figure 5. The time dependent single particle density of a harmonically trapped 
Tonks-Girardcau gas following the sudden removal of a 5 barrier of strength 
K = 10, located at position d = —6. 

the Fricdcl oscillations. At t = it/ 2 there is complete relaxation of the oscillations, 
this is followed by a complete revival at t = 7T. This is precisely the fine structure we 
see in between revivals of the density dip in Fig. 

4-2. Energy profile 

It is also interesting to consider the energy of the system as a function of the position 
of the impurity. Due to the Fermi-Bose duality, this can simply be calculated adding 
the cigcnencrgics of the single particle states up to the Fermi energy 

N-l 

E TG = ^ e« ■ ( 18 ) 

n=0 

The energy profile for a gas of thirty particles for different 5-pcrturbation strengths as 
a function of the position of the perturbation is shown in Fig. [7] One can see a series of 
exactly fifteen lobes (the plot is symmetric for d < 0.), which become more pronounced 
as the strength of the perturbation is increased, but whose position is independent of 
this increase. The position of these local maxima correspond to the positions in which 
the probabilities for the single particle wavefunctions peak, highlighting the crystal 
structure of the ground state. The increase in energy due to the impurity is largest 
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Figure 6. Detailed view of the time evolution of the single particle density for 
half a trap period. 

when it is located at the centre of the cloud. In fact, one can view the perturbation as 
a probe which, when dragged adiabatically through a Tonks gas, makes it possible to 
gain information on the crystal structure of the state through an energy measurement. 
In the next section we will look at non-local properties of the ground state properties, 
which can be calculated from the reduced single particle density matrix. 

4-3. Reduced single particle density matrix 
The calculation of the reduced density matrix, 

/oo 
f* B (x,x 2 ,.-.,x N ,t) X 
-00 

^b(x',X2, • • • , XN,t)dx 2 ■ ■ ■ dx N , (19) 

and related observables for an ultracold gas is, in general, a difficult feat. The 
dimension of the integral in eq. (|19[) increases with particle number and this is 
very demanding on computer memory resources. For a Tonks gas in a harmonic 
potential, studies have therefore mainly used numerical methods such as Monte-Carlo 
integration to calculate the RSPDM , but some analytic work has also been done 
in this direction [57]. Recently, an exceptionally efficient algorithm for calculating 
the RSPDM of a Tonks-Girardcau gas in an arbitrary external potential has been 
presented by Pezer and Buljan [28] . This algorithm allows for a numerically exact 
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Figure 7. The energy of a harmonically trapped Tonks-Girardeau gas of N=30 
particles as a function of the position of the (5-perturbation. 



calculation of the RSPDM for a large number of particles with limited memory 
resources and at a rapid computational speed. The algorithm works for both time 
dependent and time independent potentials. The essential idea is that p 1 (x,x',t) can 
be expressed in terms of the dynamically evolving single particle energy eigenbasis, 
ipi(x,t), as 



N 



p 1 (x,x',t) =2_\ij)*(x,t)Ai j {x,x' ',t)i/>j(x' ,t) . 



(20) 



The N x N matrix, A(x, x',t), is given by 

A(x,x',t) = (P -1 ) T detP , 
where the entries of the matrix P are computed as 



(21) 



(22) 



and where Sij is the Kronecker delta. Given a pair of points (a;, x 1 ) and the single 
particle basis functions xpi(x,t) one can calculate the RSPDM of a Tonks gas by 
merely calculating an N x N matrix P, its inverse and its determinant, which is a 
significant saving on computational resources. 
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Figure 8. The largest eigenvalue of the RSPDM Ao and the peak of the 
momentum distribution n(k = 0) as a function of eccentricity d for perturbation 
strengths k = 1, 5, 10 for a Tonks-Girardcau gas of N = 10 particles. 




d d 

Figure 9. The largest eigenvalue of the RSPDM Ao and the peak of the 
momentum distribution n(k = 0) as a function of eccentricity d for perturbation 
strengths k = 1,5, 10 for a Tonks-Girardcau gas of N = 11 particles. 

The RSPDM is a hcrmitian matrix and one can therefore write its spectral 
decomposition as 

p 1 (x,x',t) = Mt)#(z, • (23) 

i=0 

The eigenvectors, fa, are often referred to as natural orbitals in many-body physics and 
they provide an alternative basis of orthonormal single particle states for the many- 
body system. The eigenvalues, Xi, are known as the occupation numbers obeying the 
normalisation condition J^,- = N. The expansion in eq. (|23[) is extremely useful 
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for understanding the ground state properties of cold atomic gases, as the natural 
orbitals arc defined not only for an ideal gas but also for interacting, thermal and 
non-uniform gases. In order to obtain the eigenvectors and eigenvalues of the reduced 
density matrix, one must solve the integral value equation, 

J dx'p 1 (x,x',t)4> i (x',t) = \i(t)<t>i(x,t) . (24) 

The fraction of particles that are in the lowest lying orbital 4>o(x,t) is related to 
the largest eigenvalue Ao of the RSPDM by / = Therefore, in analogy to the 
macroscopic occupation of a single eigenstate in a Bosc-Einstcin condensate, this 
orbital is sometimes referred to as the BEC state and the quantity Ao can act as 
a measure of the coherence in the system. 

Another important measure of coherence can be derived from the momentum 
distribution of the gas which is defined as 

n(k,t) = {2-kK)- 1 J J p 1 (x,a;',t)e lfe<X R~ X l dx dx 1 (25) 
= J>(*)M 2 > (26) 

i 

where //, = {2-kK)~^ J dx e 1 ^ cf>i(x,t) are the Fourier transforms of the natural 
orbitals. For a homogenous, non-interacting Bose gas at zero temperature the 
momentum function is a (^-function, reflecting the macroscopic occupation of the 
lowest natural orbital, whereas in the strongly interacting Tonks- Girardeau gas in 
equillibrium, the momentum distribution is comprised of a non-trivial distribution of 
quasi-momenta. The amplitude of the peak of the momentum distribution at k = 
can therefore also be used to measure the spatial coherence present in the system. It 
is well known that this quantity does not follow a trivial behaviour in a disturbed, 
strongly interacting gas and that it in particular can show a dependence on having an 
even or an odd number of particles P~2| HH 1301 OH] ■ In Fig. [5] we show both quantities 
defined above, Ao and n(k = 0), as a function of distance of the perturbation from the 
trap centre for gas of TV = 10 particles. As expected, both quantities exhibit similar 
features dominated by an oscillatory structure which becomes more pronounced as 
the strength of the (5-perturbation increases. As in the interpretation of the energy 
oscillations, the maxima correspond to positions where the single particle eigenstates 
have large probabilities, i.e. where the position of the (5-potential corresponds to a 
lattice point of the underlying crystal structure of the Tonks gas. For comparison 
in Fig. [5]we show the same quantities for the N = 11 case. In this case we see a 
maximum of coherence at d = 0, as symmetry reasons require a single particle to sit in 
the centre of the trap for odd particle numbers [13] . This dependence of the coherence 
on the position of the disturbance and can be experimentally observed by measuring 
the visibility in an interference experiment. We show in Fig. [TU| the single particle 
density for a gas of N = 10 particles for which the disturbance is either located at a 
maximum of coherence (left panel), or at a minimum of the coherence (right panel). 
The difference in the visibility of the interference fringes is clearly observable. 

5. Conclusions 

In conclusion we have introduced a model to describe the harmonically trapped Tonks- 
Girardeau gas in the presence of a (5-perturbation of arbitary strength and eccentricity. 




Figure 10. Single particle density of the free evolution of the gas after realease 
from the trapping potential. Two cases are shown, (a) the delta barrier is 
positioned at d = —0.5 corresponding to a maximum coherence and (b) where 
the delta barrier is positioned at d = — 1 corresponding to minimum coherence. 
Interference fringes indicate greater spatial coherence which confirms predictions 
of Fig® 



Wc have clearly outlined the solution of the fundamental single particle problem and 
shown how it can be used in combination with known techniques to describe the 
groundstate properties of the gas. We have calculated the energy profile of the gas in 
the presence of the impurity and found an undulating profile as the perturbation 
is displaced through the gas, which highlights the crystal structure of the Tonks 
groundstate. In addition, wc have calculated the momentum distribution and largest 
eigenvalue of the RSPDM as a function of both the eccentricity and strength of the 
perturbation. Again, we have found that these properties reflect the highly localised 
nature of the particles in the Tonks gas. Furthermore as a novel application of our 
model we have investigated the time density dynamics after a sudden removal of the 
perturbation. Wc find that the gas exhibits the classical Talbot effect with the image 
of the impurity reappearing at multiples of the inverse trap frequency. Given the 
recent experiments on out-of equillibrium dynamics of a Tonks gas [9j [10], it is an 
interesting extension of our work to study the dynamics of a moving impurity. This 
is currently work in progress. 
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